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We demonstrate that perturbative QCD allows one to calculate the absolute cross 
section of diffractive exclusive production of photons at large Q 2 at HERA, while 
the aligned jet model allows one to estimate the cross section for intermediate Q 2 ~ 
2GeV 2 . Furthermore, we find that the imaginary part of the amplitude for the 
production of real photons is larger than the imaginary part of the corresponding 
DIS amplitude by a factor of about 2, leading to predictions of a significant counting 
rate for the current generation of experiments at HERA. We also find a large 
azimuthal angle asymmetry in ep scattering for HERA kinematics which allows one 
to directly measure the real part of the DVCS amplitude and hence the nondiagonal 
parton distributions. 



Recent data from HERA has spurred great interest in exclusive or diffrac- 
tive direct production of photons in e — p scattering (DVCS- deeply virtual 
Compton scattering) as another source to obtain more information about the 
gluon distribution inside the proton for nonforward scattering. In recent years 
studies of diffractive vector meson production and deeply virtual Compton 
scattering has greatly increased our theoretical understanding about the gluon 
distribution in nonforward kinematics and how it compares_to the gluon distri- 
bution in the forward direction, see list of references in Ref.EJ. We are interested 
in the production of a real photon compared to the inclusive DIS cross section. 
The exclusive process is nonforward in its nature, since the photon initiating 
the process is virtual (q 2 < 0) and the final state photon is real, forcing a small 
but finite momentum transfer to the target proton i.e forcing a nonforward 
kinematic situation as we would like. 

Similar to the case of deep inelastic scattering, in real photon production 
it is possible to calculate within perturbative QCD the Q 2 evolution of the 
amplitude but not its value at the normalization point at Qq ~ few GeV 2 
where it is given by nonperturbative effects. Hence we start by discussing 
expectations for this region. It was demonstrated inn that the aligned jet model 
coupled with the idea of color screening provides a reasonable semiquantitative 
description of F 2 n{x < 10~ 2 ,Qq). One can write a to t(l*N) using the Gribov 
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dispersion representation^ as El: 

. g f°° vtot( u AJM» - N)R e+e ~ {M 2 )M 2 ^ 
a tot {l N) = - w - M2)2 dM , (1) 

Based on the logic of the local quark-hadron duality (see e.g§ and references 
therein) we take the lower limit of integration M§ ~ Q.bGeV 2 < m 2 . In the 
case of real photon production the imaginary part of the amplitude for t = 
is obtained from Eq. |l| by replacing one of the propagators by 1/M 2 . Using 
Eq. |l| we find 

_ ImA( 7 * + N -» 7* + AQ^q Q 2 _ v a . 2 

* = i m A(r + N^ 1 + N)^ ^WTlq ln ( 1 + ^/-o)- (2) 

For A/ 2 - 0.4 0.6 GeV 2 and Q 2 as 2 - 3 GeV 2 Eq.| leads to R w 0.5. A 
similar value of R has been found in El 

The process of exclusive direct production of photons in first nontrivial 
order can be calculated as the sum of hard contributions within the framework 
of QCD evolution equationsu and a soft contribution which we evaluated above 
within the aligned jet model. In order to calculate the imaginary part of the 
amplitude, we need to calculate the hard amplitude as well as the gluon-nucleon 
scattering plus the soft aligned jet model contribution. W^find the following 
general expression for the imaginary part of the amplitude 

rQ 2 jo' 2 r 1 dr, 

ImA(x,Q 2 ,Q 2 )=ImA(x,Q 2 ) + 4Tr 2 a 2* / =±[P qg (x/x u A/xi) 

JQI Q Jx x l 

g{xi,x 2 ,Q' 2 ) + P qq (x/x 1 ,A/x 1 )q(xx,X2,Q' 2 )], (3) 

where P qg and P qq are the evolution kernels and are taken from Ref. 0. For the 
non-perturbative input, ImA(x, Q 2 ,) we will be able to use the previous aligned 
jet model analysis. We now calculate the total imaginary part of the amplitude 
from Eq. || where we used CTEQ3L for the quark and gluon distributions. 
We neglect the x 2 dependence for the moment □. The result for the ratio 
R = ImA(-y* +p — > 7* + p)/ ImAlj* +p — > 7 +p) is given in the x range from 
10~ 4 to 10~ 2 and for a Q 2 of 3.5, 12 and 45 GeV 2 relevant at HERA. According 
to our previous discussion we chose the initial distribution for the imaginary 
part of the DVCS amplitude to be twice that of the initial distribution for the 
imaginary part of the DIS amplitude. We find R to be between 0.551, 0.573 
and 0.57 for x = 10~ 4 , 0.541, 0.562 and 0.557 for x = 10~ 3 and 0.518, 0.519 
and 0.505 for x = 10~ 2 in the given Q 2 range. The ratio R will approach 
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1/2 as Q 2 is decreased to the nonperturbative scale since this is our aligned 
jet model estimate The reason for the deviation from R = 1/2 is due to the 
difference in the evolution kernels. 

As far as the complete amplitude at small x is concerned, we can recon- 
struct the real part via dispersion relations HQ which to a very good approxi- 
mation is given by: 

Re A 7r d Ya.(xImA) 
V = 7^A = 2 dlni ■ () 

X 

Therefore, our claims for the imaginary part of the amplitude also hold for the 
whole amplitude at small x. 

To check the feasibility of measuring a DVCS signal against the DIS back- 
ground, we will be interested in the fractional number of DIS events to DVCS 
events at HERA given by: 



cr(7*+p^7+p) (fcr(7* +p -> 7 + p) 1 
R ~< = aUVP) dt ' t= ° X B 1 ^ P) 

where 77 = Re A/Im A ~ 0.09 — 0.27 for the given Q 2 range, is given by 
Eq. H Note that only da/dt(t = 0) is calculable in QCD. The t dependence 
is taken from the data fits to hard diffractive processes. We computed R 7 
for x between 10 -4 and 10" 2 and for a Q 2 of 2, 3.5, 12 and 45-GeV 2 with 
the following results, where the numbers for F2 were taken from El We find 
R 1 ~ 1.1 x 10~ 3 , 4.9 x 10~ 4 at x = 10~ 4 , 10~ 3 and Q 2 = 2GeV 2 ; R 1 ~ 
1.07 x 10~ 3 , 9.3 x 10" 4 at x = 10" 4 , 10" 3 and Q 2 = 3.5GeV 2 ; R 1 ~ 4.5 x 
10~ 4 , 3.78 x 10~ 4 2.5 x 10" 4 at x = 10~ 4 , 10~ 3 , 10~ 2 and Q 2 = 12GeV 2 ; and 
finally i? 7 ~ 1.49 x 10~ 4 , 1.04 x 10~ 4 at x = 10~ 3 , 10~ 2 and Q 2 = 45GeV 2 . 

The complete cross section of exclusive photon production includes the 
Bethe-Heitler and DVCS process plus an interference term stemming from 
both processes. We find the cross section for DVCS, Bethe-Heitler and the 
interference term to be 

^ na3s Al + (l-y?)e-^Fi(x,QZ)(l + r,*), (6) 
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do bh 8a 3 sy 4 



dxdyd\t\d(f> r ir\t\Q 4 (1 - y) 



G E (t) + ^G 2 M (t) 

1 + A- 



cos 2 ( 



(7) 



d^vcs+BH _ ±2r a 3 S y 2 F 2 (x : Q 2 )co S {^)e- B \ t \l 2 hjl + jl-yf) 
dxdyd\t\d4r V RQ 5 V — 2/) 
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Figure 1: a) D is plotted versus -t for x = 10~ 4 , 1CT 3 , 10~ 2 , Q 2 = 12 GeV 2 , B = 5 GeV -2 
and y = 0.4. The solid curve is for x = 10 — 4 , the dotted one for x = 10 — 3 and the dashed one 
for x = 10~ 2 . b) The asymmetry A is plotted versus —t for x = 10 -4 (solid curve), x = 10 -3 
(dotted curve) and x = 10 -2 (dashed curve) again for Q 2 = 12 GeV 2 , B = 5 GeV" 2 and 
y = 0.4. 



G E (t) 



A-G M {t) 



1 



1*1 

l/ii-!. 



(8) 



with y = 1 — E'/E, where £7' is the energy of the electron in the final state, 
4>r = 4>n — 4>e, where 4>n is the azimuthal angle of the final state proton 
with respect to the reaction axis and </> e is the azimuthal angle of the final 
state electron and where Gs(i) and Gm(£) are the electric and nucleon form 
factors. We describe them using the well known dipole fit. The + sign in the 
interference formula corresponds to an electron scattering off a proton and the 
- sign corresponds to the positron. The total cross section is then just the sum 
of the three terms in Eq. ||[ 

At this point it is important to determine how large the Bethe-Heitler 
background is as compared to DVCS for HERA kinematics, hence, in the fol- 
lowing discussion, we will estimate the ratio D =< daovcs+BH /dxdydt > / < 
do bh I dxdydt > — 1 allowing a background comparison with <...>= L * d(f> r - 
We find that D > 1 (see Fig. |a)) for x = 1(T 4 , Q 2 = 12 GeV 2 , y = 0.4 and 
—t ~ 0.25 GeV -2 . It is convenient to illustrate the magnitude of the interefer- 
ence term in the total cross section by considering the asymmetry for proton 
and cither electron or positron to be in the same and opposite hemispheres 
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br-doDVCS+BH 



r3ir/2 
Jtt/2 ' 



b r doDvcs+BH}/ Jq 2?T d(t> r da DV cs+BH 
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Fig. [l]b shows A for the same kinematics as above and we find that the asym- 
metry is fairly sizeable already for small t . Due to this fairly large asymmetry, 
one has a first chance to access nondiagonal parton distributions through it. 
Also, we find that A is very sensitive to the energy dependence of the scattering 
amplitude since it is proportional to rj. As a result measurements of A would 
provide a window to the energy behaviour of F2n{x,Q 2 ) beyond the energy 
range currently available at HERA. 

Conclusion: In the above said we have shown that pQCD is applicable 
to exclusive photoproduction. We wrote down an equation for the imaginary 
part of the amplitude, which can be generalized to the complete amplitude 
at small x. We also found that the imaginary part of the amplitude of the 
production of a real photon is larger than the one in the case of DIS in a broad 
range of Q 2 . We also found the same to be true for the full amplitude at small 
x. We also make experimentally testable predictions for the number of real 
photon events and suggest that the number of events are small but not too 
small to be detected at HERA. Finally, we demonstrated that first the Bethe- 
Heitler process does not overshadow the DVCS cross section and secondly 
that measuring the asymmetry A at HERA would allow one to determine the 
real part of the DVCS amplitude, i.e. , gain a first experimental insight into 
nondiagonal parton distributions. 
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This effect will be taken into account in the actual numerical calculation 
by incorporating the different evolution of diagonal and nondiagonal dis- 
tribution in the formula for the imaginary part of the amplitude. 
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